Estimation of steady-state basic parameters of stars 
B.V.Vasiliev 

Institute in Physical and Technical Problems, Dubna, Russia, 141980 

O 

q . e — mail : vasiliev@mainl.jinr.ru 

(N 

a. 

Abstract 

' From a minimum of total energy of celestial bodies, their basic 

parameters are obtained. The steady-state values of mass, radius, and 
, temperature of stars and white dwarfs, as well as masses of pulsars are 

calculated. The luminosity and giromagnetic ratio of celestial bodies 
are estimated. All the obtained values are in a satisfactory agreement 
with observation data. 
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j^' 1 Introduction 

o 

It is agreed-upon to consider that pressure and density of a substance inside 
stars grow more or less monotonously approaching their centrs. It is supposed 
even that an ultrahigh-density substance can exist near the center of some 
^ ■ stars. Thus, it is considered that the basic parameters of stars, such as mass, 

radius and temperature, are more or less arbitrary. 

However, density, mass, and temperature of a star can be determined 
coming from the requirement of a minimum of energy if a star exists in an 
equilibrium state. 

The effect of gravity- induced electric polarization of plasma inside a star 
gives a basis for such a calculation. 

Under the influence of gravity, the nuclei of plasma "hang" on their elec- 
tronic clouds. Therefore, they are always located right below these clouds. 
As a result, electrons form a stratum at the external surface of a plasma 
body, and inside a star each cell of plasma obtains a positive charge 8q. This 
charge is very small and equal to jp-e ~ 10~ 19 e by the order of magnitude (oq 



1 



is Bohr radius, R s t is the radius of a star, and e is the charge of an electron). 
However, it is sufficient, that at each point inside a star a gravitational force 
is cancelled by an electric force 

mi g + 5qE = } (1) 

where rrii is the mass of an ion, g is gravity acceleration, E is an electric field 
intensity created by charges Sq 

divE = AirSqn, (2) 

where n is the electron density. 

From here it follows that the gradient of pressure inside a star is equal to 
zero, and the substance has a steady value of density. For a star consisting of 
degenerate relativistic plasma, the temperature effects can be neglected. For 
a star consisting of hot nonrelativistic plasma, the absence of the gradient of 
pressure and density of a substance requires an assumption that temperature 
inside a star is constant. It concerns, certainly, only the basic core of stars 
consisting of electrically polarized plasma. There is a thin stratum at the 
surface of a star, where a substance exists in an usual atomic state. In this 
stratum, the gradient of pressure, the gradient of density, and the gradient of 
temperature are present. Because of a small thickness of this stratum (about 
10~ 2 R st ), its role in the energy balance of the whole star can be neglected. 

Since the gravitational force is cancelled by the electric force in cores of 
stars, the gravitational energy is compensated by the electric energy inside a 
core as well |IJ (the gravitational energy is negative). There is an uncompen- 
sated part of energy of a gravitational field. It is the energy of gravitational 
field outside a star 

Eg - -- ^ , (3) 

where G is the gravitational constant, M and R are mass and radius of a star. 
Because plasma is electroneutral as a whole, one proton should be related 
to electron of the Fermi gas of plasma. The existence of one neutron per 
proton is characteristic of a substance consisting of light nuclei. The quantity 
of neutrons increases approximately to 1.8 per proton for the heavy- nuclei 

substance, i.e. to one electron in plasma a mass equal to m p (^j should be 

related. Thus 
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M = Nm p (^y (4) 

Here N is a full number of electrons in a star, m p is the proton mass, A and 
Z are the mass number and charge of a nucleus, respectively. For hydrogen 

jj^J — 1. Therefore, 

!<(!)< 2 - 8 - (5) 

The field Eq.(|j) tends to compress a star as a whole. 

The jump in electric polarization at the surface of a star core produces 
additional pressure inside a star [U, since the jump in polarization is accom- 
panied by the pressure jump. As this energy is considerably smaller than the 
gravitational one, it can be neglected for simplicity. 

At high temperature the energy of the black radiation can be important 
in a general balance of the star energy 

4rr 

E r = —T A V (6) 

c 

where a = 5.7 • 10 -5 — fur is the Stefan-Boltzmann constant, and V is the 
volume of a star. 

Thus, the total energy of a star is a sum of the total energy of plasma 
and the energy of black radiation 

Etotai = E p i + E r . (7) 

2 A star consisting of hot nonrelativistic elec- 
tron - nuclear plasma 

The energy of plasma in a star is a sum of its potential energy U and kinetic 
energy E^ inetic . According to the virial theorem the potential energy of a 
system of interacting charged particles and their kinetic energy are connected 



U — —2Ekinetic- (8) 
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Because the ion mass is large, the chemical potential of ions is small and 
their kinetic energy can be neglected. Thus, electrons will make a main 
contribution into the kinetic energy of plasma. 

We name plasma a hot one, if the temperature inside a star is more then 
the degeneration temperature of the electron gas and the electric interaction 
between particles of plasma can be neglected 



e 2 



(9) 



Under this condition it is possible to consider the electron gas of plasma an 
ideal one, and write its equation of state as the ideal gas law 

PV = NkT. (10) 
Thus, the kinetic energy of plasma in a star 

3 

E kinetic = 2 NkT > 

where k is the Boltzmann constant. 

Because gravity does not influence photons and acts on heavy ions only, 
the total energy of hot plasma is 

E p i = _|_ -NkT, (12) 

pl 2R 2 K J 

or, according to the virial theorem (Eq.(^D), the energy of plasma in a 
star 

E " = "TiT = J 2 NkT (13) 

and the total energy of a hot star is 

GM 2 3 ArrrT1 4(7 . r GM 2 4cr . . ^ 

Etotai = + -NkT + — T 4 V = + — T 4 V 14 

total -r -r -r V I 
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3 Calculation of the corrections to the elec- 
tron energy 

This equation is valid for an electron gas at very high temperature, when the 
identity of particles is of no significance. At finite temperature the identity 
of particles leads to an extra stiffness of the electron gas than it follows 
from Eq. flTOl) . On the other hand, electron gas in plasma will exhibit more 
softness under the influence of a Coulomb field of nuclei. We shall carry 
out the account of these corrections according to the methods described by 
Landau and Lifshits 0. 

3.1 The correction on electron identity 

The total energy of the electron gas is known 



where e = P 2 /2m, P and m are momentum and mass of an electron in the 
nonrelativistic gas. 

According to the definition of a chemical potential of electrons || , 



Keeping the first terms of series, an approximated expression of the free 
energy of the hot electron gas with the account of the correction on the 
identity of electrons is obtained 




(15) 




(16) 



and at high temperature 




F = F ldeal + N 



7r 3 / 2 e 3 a^ 2 n 
4(fcT) 1 /2 



(18) 
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3.2 The estimation of the influence of nuclei on the 
hot electron gas 

The nuclei reduce pressure of electron gas apart from them, therefore, this 
correction to pressure must be negative. The estimation of this correlative 
correction can be carried out by the method offered by Debye and Hiickel 

§• 

The energy of an charged particle inside plasma is equal to e<p e , where 
e is a charge of a particle, and tp e is the electric potential created by other 
particles on this particle. This potential in plasma is determined by the 
Debye law 

m = -e r o . (19) 
r 

As in dense plasma nuclei form an ordered lattice @, electrons play a role in 
screening charges only. 

Thus, the Debye radius is 



/ kT 

At the small value of ratio the potential can be expanded into the series 



*>=--- + ... (21) 
r ro 

The following terms are converted into zero at r = 0. The first term of 
this series is the potential of the electron. The second term is the potential 
created by other particles on electron. Therefore, the energy of the electron 
gas in plasma created by the action of screened fields of other particle: 

e 2 „ .3 nriv 3 / 2 



E corr = N = 2e 3 J rn r. (22) 



By using of the known Gibbs-Helmholtz equation 

E d F 



7 ~ O'V I 



we obtain the correction to the free energy 
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3.3 The steady-state value of a star density 

Accounting for these corrections on electron identity and influence of ions, 
the free energy of the hot electron gas acquires the form 

3/2 3 3/2 9 3 1/2 

f = ^ + ^»-4^ (25) 

From the equilibrium condition 

we obtain a steady-state value of electron density in plasma of a star 

n st = 7T~2~3 - 2 ■ 10 23 cm" 3 . (27) 

For a star with a sufficiently high temperature, according to Eq.(^), it is 
energetically favourable to have this density 

3.4 The steady-state values of the mass, the radius and 
temperature of a star 

According to Eq. (^4|) , the total energy of a star 

-kTn st V+^ 
2 c 

This function has a minimum at the temperature 



toted — -~kTn st V + —T A V. (28) 



'3cfcn,A 1/3 / ck \ 1 / 3 

According to the virial theorem Eq.(||), at an equilibrium state into the 
electron gas subsystem 



GM 2 , s 

— = 3NkT. (30) 
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With the account of Eq.(^) and Eg. fl2"9|) , we can calculate a steady-state 
value of radius of hot stars 



R 



3^ 
2 
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A8aG 3 m 6 p 



1/6- 



a 



A- 
■ Z' 



1.6 

(!) 



Re. 







where R & is the radius of the Sun. 

The steady-state value of mass of hot stars 



(31) 



M, 



27 



r. 



4 V37raG 3 



3/2 



(4) 2 



where M is the mass of the Sun. 



(32) 



4 The equilibrium condition of a white dwarf 

When density of a substance is sufficiently high and temperature is insignifi- 
cant, the electrons form a degenerate and relativistic Fermi gas. Thus, a star 
obtains the other equilibrium state. 

It is known that the electron g cLS IS cL relativistic one, when its Fermi 
momentum 

p F = (3n 2 ) 1/3 n 1/3 h > mc, (33) 

i.e., at density of particles n > 10 31 cm -3 , which is characteristic for white 
dwarfs. 

For a degeneration, electron gas temperature must be 

T <<: ~T~ ~ 10l ° K - ( 34 ) 
The energy of the relativistic electron gas is known 

E^( 3 A\ cN ay<\ (35) 



4V J \V 
Therefore, the total energy of a white dwarf 



GM 2 3/ 2 \V3 /JV\V3 

E ^ = -^ + i^) hcN [v) ■ (36) 
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As temperature is defined by the inequality Eg. (|34T) , the temperature de- 
pending terms can be neglected. It is possible under the condition 

—T 4 V « - A {^ 2 ) l/ \cN( ^V /3 . (37) 



c 4 V ' \V 

As according to the virial theorem (Eq. 



gm 2 3/ 2 y/\ /iV\ 1/3 . . 

- Svr 2 hcN(-) , (38) 



2R 2 V J \V 
the full number of electrons inside any white dwarf is fixed 

M 3 5 /V/ 2 / ftc \ 3 / 2 3.2 -10 58 



Nda ~mm P - 2 U(i)^J * '(if • (39) 



Thus, in the equilibrium state any white dwarf should have the steady-state 



value of mass, depending on world constants only and the ratio ( 4 ) : 



M, 



(I w 



3 5 /V/ 2 / he \ 3 / 2 



\Gm 2 



rr^_ ^ 12.5M 

(A\2 — ( A\2 ■ 



According to the inequality Eq . (|33"D , the steady-state value of radius of a 
white dwarf 

/Q,7T \ 1/3 fj 

R dw < ~rN dw — ~ 3 • 10"=% (41) 
V 4 / mc 



and the density of electrons 



n dw = » 10 30 cm- 3 . (42) 

3 K'dw 



This correlates well with Eq. (|33"D . 
According to Eq.(^) 



mc 3/2 



T dw « — m ~ 3 ■ 10 9 K. (43) 

(4vr) 1 /2 a 1 /^ 3 / 4 v ; 



Which is consistent with Eq.( 

The steady-state value of mass of a star (|40|), consisting from the rela- 
tivistic electron gas, does not depend on mass of an electron. Due to this 
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fact, this equation is valid for stars, consisting of a degenerate gas of other 
relativistic Fermi particles, in particular, neutrons. Therefore, if we can con- 
sider a pulsar as a neutron star containing a small number of protons and 
electrons (with concentration more than 10~ 19 ), the equation Eq.(|40"D should 
determine the steady-state value of mass of a pulsar. 



5 The mass- luminosity relationship 

It seems possible to assume that temperature inside a star determines tem- 
perature at its surface. If the relation between a surface temperature T and 
the steady-state values of inner temperature T st is constant for all stars 

T = const ■ T st , (44) 

it is possible to deduce a relationship between a visible luminosity of a star 
and its other parameters, in particular, mass. 
As the luminosity of a star 

L = AtcR 2 -T\ (45) 
n 

using the obtained above expressions and after simple transformations, we 
obtain 

L ( M \ 10 / 3 , . 

(46) 



where L & is the luminosity of the Sun. 



6 The comparison of the calculated values 
with the observation data 

The masses of stars can be measured with a considerable accuracy, if these 
stars compose a binary system. There are almost 300 double stars which 
masses are known with the required accuracy 0]. Among these stars there 
are stars of the main sequence, dwarfs, and pulsars. Approximately one half 
of them are visual binaries. Their masses are measured with the a high pre- 
cision. Other half consists of spectroscopic binaries and eclipsing binaries. 
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For these stars the accuracy of mass measurement is a slightly worse. Never- 
theless, we shell carry out a comparison of the calculated parameters of stars 
on the basis of all these measurements. 



6.1 Masses of celestial bodies 

According to these data, the distribution of masses of hot stars is described 
by the equality 

< M >= (2.98 ± O.25)M . (47) 

It is in satisfactory agreements with the calculated steady-state value of mass 
of a star Eq.(^) by the order of magnitude. Graphically this distribution is 
shown in Fig.l. According to the binary star tables 0], the distribution of 
white dwarfs masses is described by 

< M >= (0.96 ± O.O5)M . (48) 

This distribution is shown in Fig. 2. In this figure, the distribution of pulsar 
masses || is shown. It is described by 

< M >= (1.40 ± 0.02)M Q . (49) 

It is possible to conclude that with the account of the corrections on the 
factor J , the results of calculations of the steady-state values of masses 
are in agreement with measured data. 



6.2 The mass-luminosity relationship 

The luminosity of stars from binary systems, depending on their masses are 
shown in Fig. 3. The results of more precise measurements of visual binaries 
are marked by squares. The data for spectroscopic and eclipsing binaries are 
marked by triangles and points, respectively. The line corresponds to the 
calculated dependence Eq . (|46|) . It can be seen from this figure, that in the 
logarithmic scale the calculated dependence is in satisfactory agreement with 
the observed data. 
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7 The gyromagnetic ratio of stars 



The effect of gravity-induced electric polarization JT| is characteristic for all 
celestial bodies, consisting of a substance in a plasma state. The value of 
this polarization does not depend on plasma properties: no matter if it is a 
relativistic one or not and if it has a degeneration or not. As divE = Anp 
and divg = —AnGj and according to Eq. ([!]), the gravity-induced density of 
the volume electrical charge p is related to the density of substance 7 by the 
relation [[!]]: 



Therefore a rotation of a celestial body about its axis must induce a magnetic 
field. If one assumes that the electrically polarized core of a body occupies 
its entire volume and if one neglects the existence of a surface stratum, where 
the substance is in an atomic state, that the gyromagnetic ratio (the relation 
of a magnetic moment of a body to its angular momentum) will obtain a 
steady value [|TJ: 



It is possible to consider the above assumption is quite acceptable for large 
celestial bodies, such as stars, and less acceptable for planets, where the 
mantle can be rather large. However, the detailed calculation for the Earth 
gives the result which is in agreement with the measured value with in the 
accuracy of the factor of two |4|. The values of giromagnetic ratio for all 
celestial bodies (for which they are known today) are shown in Fig. 4. 

The data for planets are taken from 0, the data for stars are taken from 
H, and for pulsars - from ||. Therefore, for all celestial bodies - for planets 
and their satellites, for Ap-stars and several pulsars - the calculated value 
Eq.(|5i~D with a logarithmic precision quite satisfactorily agrees to measure- 
ments, when moments itself change within the limits of more than 20 orders. 



p = Vgj. 



(50) 




(51) 
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Figure 1: The mass distribution of stars from the binary systems ||. Data 
for visual, spectroscopic, and eclipsing binaries are shown separately. On the 
abscissa, the logarithm of the star mass over the Sun mass is shown. 
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Figure 2: The mass distribution of white dwarfs and pulsars from the binary 
systems || [|. On the abscissa, the logarithm of the star mass over the Sun 
mass is shown. 
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Figure 3: The luminosity of stars from binary systems, depending on their 
masses are shown. The results of more precise measurements of visual bina- 
ries are marked by squares. The data for spectroscopic and eclipsing binaries 
are marked by triangles and points, respectively. The line corresponds to the 
calculated dependence Eq.(|46|). On the ordinate, the logarithm of the star 
luminosity over the Sun luminosity is shown. On the abscissa, the logarithm 
of the star mass over the Sun mass is shown. 
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Figure 4: The observed values of the magnetic moments of celestial bodies 
vs. their angular momenta. On the ordinate, the logarithm of the magnetic 
moment over Gs ■ cm 3 is plotted; on the abscissa the logarithm of the angular 
momentum over erg ■ s is shown. The solid line illustrates Eq.(^l|). The dash- 
dotted line is the fitting of the observed values. 
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